Decoherence and Dynamical Entropy Generation in Quantum Field Theory 
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We formulate a novel approach to decoherence based on neglecting observationally inaccessible 
correlators. We apply our formalism to a renormalised interacting quantum field theoretical model. 
Using out-of-equilibrium field theory techniques we show that the Gaussian von Neumann entropy 
for a pure quantum state increases to the interacting thermal entropy. This quantifies decoherence 
and thus measures how classical our pure state has become. The decoherence rate is equal to the 
single particle decay rate in our model. We also compare our approach to existing approaches to 
decoherence in a simple quantum mechanical model. We show that the entropy following from the 
perturbative master equation suffers from physically unacceptable secular growth. 
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I. INTRODUCTION 

Decoherence is a quantum phenomenon that describes 
how a quantum system turns into a classical stochastic 
system. The theory of decoherence is widely used, e.g. 
in quantum computing black hole physics 0, infla- 
tionary perturbation theory Q and in elementary parti- 
cle physics such as electroweak baryogenesis models [1]. 
It is hard to study decoherence in quantum field theory 
(QFT) as it requires out-of-equilibrium, finite temper- 
ature and interacting quantum field theoretical compu- 
tations. Let us first discuss the conventional approach 
to decoherence 0] and its shortcomings. We consider a 
system S in interaction with an environment E. The 
observer's inability to detect the environmental degrees 
of freedom allows us to trace over E to construct a re- 
duced density operator from the full density operator 
Prcd = Trg [p] . The unitary von Neumann equation for the 
density operator idtp= \H,p] is replaced by a perturbative 
"master equation" for p^cd that is no longer unitary: 



idtPrcd = [Hs,Prcd] +T^[Prcd] , 



(1) 



where Pf/Oiod] collectively refers to all non-unitary dis- 
sipative terms [1, 0] ■ If Pred is approximately diagonal 
in the pointer basis, our quantum system effectively be- 
comes a classical stochastic system. In reality, however, 
systems do not fully decohere so in order to quantify the 
amount of decoherence, one is interested in the entropy 
increase due to the loss of quantum coherence. 

The master equation ([1]) suffers from both theoretical 
and practical shortcomings. It is unsatisfactory that p^-cd 
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evolves non-unitarily while the underlying QFT is uni- 
tary. On the practical side, the master equation is so 
complex that basic field theoretical questions have so far 
never been properly addressed: no well-established treat- 
ment to take perturbative corrections to pred into account 
exists (for cases in quantum mechanics (QM) see however 
Q), nor has any piod ever been renormalised. Moreover, 
as we will show, the perturbative master equation leads 
even in very simple QM situations to physically unaccept- 
able secular growth of the entropy, which is caused by the 
perturbative approximations used in specifying ([T]). 

The goal of this letter is twofold. Firstly, we present 
our novel approach to decoherence that does not suffer 
from the shortcomings mentioned above and apply it to 
a very simple QM toy model of A^-|-l coupled simple 
harmonic oscillators (SHO's). Secondly, wc apply our 
approach to an interacting QFT, using non-equilibrium 
field theory techniques to calculate the time evolution 
of the Gaussian von Neumann entropy. The focus and 
novelty of this letter is not on calculational results, but 
instead on a clear and succinct presentation of our ap- 
proach to decoherence and its relation to entropy genera- 
tion in relativistic quantum field theories. The results of 
various calculations we discuss here serve mainly to test 
and illustrate our approach to decoherence. We refer the 
reader for further details to 



II. ENTROPY AND CORRELATORS 

It is of course widely appreciated that entropy can be 
generated as a result of an incomplete knowledge of a sys- 
tem. The infinite hierarchy of irreducible n-point func- 
tions completely captures all properties of a system, how- 
ever only a finite subset can be probed experimentally. 
This leads us to propose [^-[ll|: neglecting observation- 
ally inaccessible correlators will give rise to an increase 
in entropy of the system as perceived by the observer. 
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So both S and E evolve unitarily, however the observer 
would say that S evolves to a mixed state with positive 
entropy as information about S is dispersed in inacces- 



sible correlation functions. Giraud and Serreau |l3j ad- 
vocate similar ideas which they illustrate by an analysis 
in A0^, building on earlier work 0, The total von 
Neumann entropy can be subdivided as follows: 



5, 



vN 



5-8 (t) + 5ng(t) = sl + S% + SI;e + S'^^ . (2) 



Here, is the total Gaussian von Neumann entropy, 
that contains information about both S and E and the 
their correlations SE, and S*"^ is the total non-Gaussian 
von Neumann entropy (consisting of S, E and SE contri- 
butions). In the following, we assume that the relevant 
properties of quantum systems are encoded in the Gaus- 
sian part pg of the full density operator p. This is justified 
as higher order non-Gaussian n-point functions are per- 
turbatively suppressed. Our approach can be improved 
if e.g. three- or four-point correlators are experimentally 
accessible such that knowledge of these correlators can 
be included in the definition of the entropy [l^l- Al- 
though 5'vN is conserved in unitary theories, Sg(t) can 
increase at the expense of other decreasing contributions 
to the total von Neumann entropy, such as S'"^. Calzetta 
and Hu [l5| prove an iJ-theorcm for a quantum mechan- 
ical 0(7V)-model and refer to "correlation entropy" what 
we would call "Gaussian von Neumann entropy" . For 
e.g. the Gaussian system density matrix for a real scalar 
field, one can straightforwardly find its associated von 
Neumann entropy [TO, [13 ■ 

S'| = -Tr[/5glnpg]: 

^ ^ (3) 

Here, A is the phase space area that the sytem's state 
occupies in Wigner space {fi=l, k=\\k\\): 
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A2(fc,t)=4[(02)(fr2)_({0,7r}/2)2 



(4) 



The most efficient way of solving for the three Gaussian 
correlators in (|H) is to solve for the statistical propagator 
F^{x;y)='TT[p{to){<j){x),(t){y)}/2], from which the corre- 
lators can be straightforwardly extracted: 

{${x,t)${ij,t)) = F^{x,t;y,t%^t' (5a) 

{T:{x,t)n{y,t)) = dtdt'F^{x,t;y,t')\t=t' (5b) 

{{^{x,t),7r{y,t)}/2) ^ dt'F^{x,t;y,t')\t^t' ■ (5c) 

For a pure state we have A = l and 5| = 0, whereas for 
a mixed state A > 1 and 5*1 > 0. The phase space area is 
conserved by the evolution d/dt[A^]=0 in free theories. 
As Sg is the only invariant measure of the entropy of a 
Gaussian state [1^, we argue that 5*1 should be taken 
as the quantitative measure for decoherence. This is to 
be contrasted with most of the literature where different, 
non- invariant measures are proposed • 

Finally, let us connect the standard approach to ours. 
The expectation values in ([5]) can in principle be cal- 
culated from the reduced density matrix, consider e.g. 
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Figure 1: Exact Gaussian von Neumann entropy Sg (solid) 
and the entropy from the master equation ^ (dashed). We 
use cj„/a;o = l + n/100, l<n<50, A/w,^ = 3/40 and /3ajo = 2. 



(0^) = Trs_(_£;[/5(^^] =Tr5[/5rcd0^]- Suppose one has been 
able to solve for p in full generality, one would then be in- 
terested in evaluating a "reduced von Neumann entropy" : 
Sl''^{t)^~TT[p,^d{t)logPrcd{t)]. We thus conclude that 
in the Gaussian approach we also advocate, one would 
find Sl'^ = Sg. The main point of this letter is that solv- 
ing for a renormalised Erf, in an interacting QFT is more 
tractable than solving for a full p, as we have recently wit- 
nessed a further development of sophisticated techniques 
in out-of-equilibrium QFT [l7| . 



III. QUANTUM MECHANICS 

Let us consider the well-known SHO model of iV- 
coupled oscillators: 



L = -x-' 
2 



1 2 2 
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1 2 2 



), (6) 



where x plays the role of S, and the {qn} of E. This 
model provides us with a level playing ground to compare 
the two approaches to decoherence as QM toy models like 
these are free from the drawbacks previously mentioned: 
we can actually solve ([1]), and we do not have to worry 
about renormalisation. Let us start by considering our 
approach. Equation ([2]) simplifies to SvN = Sg + S^ + Sg^ 
as 5"s = 0. The role of the neglected non-Gaussianities 
in the QFT case is played in the quadratic QM model 
by the correlation entropy 51^. Let us consider A^=:50 
oscillators, where all E oscillators with frequencies {a;„} 
act as a thermal bath at temperature A:b = 1- Ini- 
tially, we require that 5 is in a pure state 51 = 0. Fol- 
lowing our approach to decoherence, we can numerically 
study the exact evolution of the three Gaussian QM ana- 
logues of equation ([S]): {x^), {pD and {{x,Px})/2. In the 
resonant regime where one or more a;„ ~a;o, the {qn} cou- 
ple effectively to x such that we expect to observe swift 
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thermalisation and decoherence. In figure [T] we show the 
resulting evolution of the entropy (solid line) . The oscil- 
lations are a manifestation of Poincare's recurrence the- 
orem The dot-dashed line is the thermal entropy 
corresponding to a temperature /Sujq = 2 indicating per- 
fect thermalisation, which at high temperatures, /3wo <d, 
would correspond to perfect decoherence. Given the weak 
coupling A /u!q =3/40 wc observe imperfect decoherence 
(incomplete thermalisation). If we increase A, we observe 



full thermalisation of S" 111. We conclude that a finite 



amount of entropy is generated at the expense of obtain- 
ing a negative contribution to the correlation entropy. 

This behaviour is to be contrasted with the dashed 
line, the master equation ([ij. The entropy resulting 
from evolving prcd suffers from physically unacceptable 
unbounded secular growth. In cases where the entropy 
obtained from the master equation settles to a constant 
value before it breaks down, it generically overshoots. 
Let us stress that we are in the deep perturbative (and 
resonant) regime and that all eigenfrequencies following 
from © are positive. The perturbative master equation 
incorrectly resums self-mass contributions, unlike the 2PI 
(two particle irreducible) scheme available in QFT. Away 
from the resonant regime, where all {w„} significantly dif- 
fer from ujQ, the evolution of the entropy from the master 
equation agrees up to the expected perturbative correc- 
tions with our approach, but this is the regime where 
only a negligible amount of entropy is generated (no de- 
coherence). We conclude that decoherence rates at early 
times can be calculated accurately using the standard ap- 
proach to decoherence however the perturbative master 
equation breaks down at late times such that the total 
entropy increase, i.e.: the total amount of decoherence, 
cannot be calculated reliably. 



IV. QUANTUM FIELD THEORY 

We will consider the following interacting scalar QFT: 

S[cl),x] = J d°x{/:o[0] + Coix] + A„t[0,x]} , (7) 
where: 

AM = -i9^0(x)a,(/.(x)r;^'^ - im2 02(2.) (ga) 



(8b) 



Here, will play the role of S and x of i?. We as- 
sume that {(j)) = = (x) , which can be realised by suit- 
ably renormalising the tadpoles. Equation ([2]) reduces 
to 5vn = 5|-|-S'|-K5"« as 5*1^ =0. The analogy with the 
simple QM model is clear: in the QM case 5*1 increases as 
information is neglected in the Gaussian SE correlations, 
whereas in the QFT model Sg increases as non-Gaussian 
information is neglected, mainly in the SE correlations. 
Note that in 
glected. 



131 no n- Gaussian S correlations are ne- 
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Figure 2: as a function of time (solid), which quickly settles 
to S'nis (dashed), where "ms" is an abbreviation for mixed 
state. We use I3m^ — 1/2, k/m^ = l, h/m^ = 3 and A'' = 2000 
up to tm^ = 100. 



Wc work in the Schwingcr-Kcldysh or in-in formalism. 
The 2PI effective action captures the effect of perturba- 
tive loop corrections to the various propagators iA°^ and 
iA°^, {a,b} — ±. Varying the 2PI effective action yields 
the Kadanoff-Baym (KB) equations for the Wightman 
propagators: 

{dl-ml)2Af{x;y) -Y^cf d°iMf (:r ;i)zA;=F(i;y) = 0, 

(9) 

where the self-masses for 4> at one loop read: 



iMl\x-y) 



>Af(^;y))' 



(10) 



Using dimensional regularisation, we find that a local 
mass counterterm rcnormalises our theory Q: 



z/i2r(f - i)/i°-'' 

16tt^ {D ~ 2,){D - A) 



i5"{x^v). (11) 



The 
written 



KB 
in 



equations can be solved when 
terms of the causal propagator 
jA^(a;;y) = Tr(p(to)[0(x),0(y)]) = zA^+-iA+- and 
the statistical propagator = (iA^^ + «A^^)/2. Ini- 
tially, at i = to, is in a pure state and % is in thermal 
equilibrium at temperature We assume that \'3>h 

such that we can neglect the backreaction from S on 
E. Perturbatively, this argument is justified as the 



backreaction contributes only at 0{h'^/ujj^ 



Here, a single (double) solid black line corresponds to a 
free (resummed) (f> propagator. Appreciate that the lead- 
ing order self- masses yield Gaussian corrections to lA'^'', 
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Figure 3: The interacting thermal entropy Sms as a function 
of h/m^. The dashed black line indicates the free thermal 
entropy. We use Pm,^ = l/2 and k/m^ = l. 



SO all non-Gaussianities are indeed perturbatively sup- 
pressed. In our approximation scheme x influences 4> thus 
via the IPI self-masses which are calculated elsewhere 
• Memory effects play an important role in non- 
equilibrium QFT and stem from the time integrals J_^dt 
over the self-masses in © with some cancellations such 
that the final result is causal. We deal with the initial 
time divergences by including memory effects at times 
t<to, and an appropriate truncation method of these in- 
tegrals (for alternative approaches see [l^, e.g. Garny 
and Miiller impose non-Gaussian initial conditions). 

The resulting evolution for the entropy is shown in 
figure [2] which follows from Fcj, and equations ([SHSl)- Ini- 
tially, we have Sg (tp ) = which then rapidly increases 
and settles to its asymptotic value Sms- We conclude 
that although a pure state remains pure under unitary 
evolution, the observer perceives this state over time 
as a mixed state with positive entropy as non- 
Gaussianities are dynamically generated. In other words, 
a realistic observer cannot probe all information about S 
and thus discerns a loss of coherence of our pure state. 

We can extract two relevant quantities: the maximal 
amount of decoherence and the decoherence rate. The 
maximal amount of decoherence Sms can be calculated 
by independent means and offers a powerful check of our 
approximation scheme [1, At late times, is time 
translationary invariant such that the entropy is time 
independent, which can be appreciated from equation 
e.g. F^{k,0) = /dfc7(27r)F0(fc^) is time independent. 
F^{k^) follows directly from the KB equations and the 
self- masses in Fourier space [1, [l^- We show in figure 
[3] that 5'nis reduces to the (free) thermal entropy when 
/i-^O given by Afroc = coth(/3a;0/2) and equation We 
thus see that Sms is the interacting thermal entropy. 



Figure 4: Approach to Ams (solid) and the decoherence rate 
Fdec (dashed). We use the same parameters as in figure [2] 

The decoherence rate is the rate at which the phase 
space area (or entropy) changes. We observe as leading 
order effect an exponential approach to Ams in figure H) 

SAkit) = Ams - Afe(t) ~ {Ams - Afc(O)) exp[-rdcci] • 

Here, Tdoc is the decoherence rate of our model which 
is well described by the single particle decay rate of our 
interaction: T^^^^ = — Im(zAfp/a;0, where iM^ is the 
retarded self-mass projected on the quasi particle shell: 



dec - 



32ttuj^ 16nk/3uj^ 
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V. CONCLUSION 

We have studied the time evolution of the Gaussian von 
Neumann entropy in an interacting, out-of-equilibrium, 
finite temperature QFT in a renormalised and perturba- 
tive 2PI scheme. We have extracted two relevant quanti- 
tative measures of decoherence: the maximal amount of 
decoherence and the decoherence rate. This study builds 
the QFT framework for other decoherence studies in rel- 
evant situations where different types of fields and in- 
teractions are involved. In cosmology for example, the 
decoherence of scalar gravitational perturbations can be 
induced by e.g. fiuctuating tensor modes (gravitons) 0], 
isocurvature modes [l^ or even gauge fields. In quantum 
information physics it is very likely that future quantum 
computers will involve coherent light beams that inter- 
act with other parts of the quantum computer as well as 
with an environment For a complete understanding 
of decoherence in such complex systems it is clear that a 
QFT framework such as developed here is necessary. 
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